Quantum Hall activation gaps in bilayer graphene 
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Abstract 

We have measured the quantum Hall activation gaps in bilayer graphene at filling factors v = ±4 and v = +8 in high 
magnetic fields up to 30 T. We find that energy levels can be described by a 4-band relativistic hyperbolic dispersion. 
The Landau level width is found to contain a field independent background due to an intrinsic level broadening and a 
component which increases linearly with magnetic field. 
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1. Introduction 

Bilayer graphene consists of two coupled one atom 
thick carbon layers arranged in a honeycomb lattice. 
These two layers are coupled in an A1-B2 stacking, 
where A and B represent the two triangular sublattices 
of the individual graphene layers fljj, [2J], see inset in 
Fig. Q] 

Bilayer graphene attracted the attention of re- 
searchers in transistor electronics due to its tunable gap 
induced by an electric field 121 |4j, |5[] or by doping |g]. 
At the same time its hyperbolic dispersion of massive 
chiral Dirac fermions J7[] makes it an intriguing ob- 
ject for more fundamental physics as a connecting link 
between the 2D single layer graphene with relativistic 
massless charge carriers and the 3D classical Fermi gas 
in graphite. 

In particular, placing bilayer graphene in a mag- 
netic field leads to an unconventional integer quantum 
Hall effect |2j, |8[] distinctively different from the rela- 
tivistic half-integer quantum Hall effect in single layer 
graphene Qii and the integer quantum Hall effect in 
the conventional 2D electron systems [11]. 

In this Communication we probe the hyperbolic dis- 
persion of bilayer graphene by measuring the activation 
gaps between neighboring Landau levels. We will show 
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Figure 1: (color online). Hall resistivity p„, and longitudinal resistiv- 
ity p xx as a function of magnetic field at T = 24 K for a hole-doped 
device (n = -3.07 X 10 16 m~ 2 ). The inset shows the A1-B2 stacking 
in bilayer graphene. 



that the activation energies are described by the Lan- 
dau level spectrum for a 4-band relativistic hyperbolic 
dispersion with a Landau level width composed of two 
components, namely a field independent intrinsic broad- 
ening and a linearly field dependent broadening. 



2. Experimental results and discussion 

We have measured the resistivity of a 1 /im wide Hall 
bar of bilayer graphene as a function of magnetic field 
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up to 30 T. The sample was fabricated by microme- 
chanical exfoliation of a flake from natural graphite, de- 
posited on a 300 nm SIMOX wafer and identified us- 
ing an optical microscope 11211 . The carrier concentra- 
tion in such samples can be changed by applying volt- 
age to the underlying conducting silicon, which acts as 
a back gate. To remove surface impurities we have an- 
nealed the sample at 390 K in vacuum prior to mea- 
surements, giving on an increasing the mobility from 
fx = 1000cm 2 /(Vs)tOjU = 3000 cm 2 /(Vs). 

Figure Q] shows typical traces of the Hall resistivity 
Pxy and the longitudinal resistivity p xx for a hole-doped 
sample (charge carrier concentration n = -3.1 x 10 16 
rrT 2 ) as a function of magnetic field at a temperature of 
24 K. The Hall resistivity exhibits plateaus at filling fac- 
tors v = -4/ where i = 1, 2, 3, . . . and the "— " sign indi- 
cates holes. These plateaus are accompanied by minima 
in the longitudinal resistivity. 

The appearance of quantum Hall plateaus at v = ±4i 
can be understood from the Landau level spectrum of 
bilayer graphene shown in Fig.|2ja). The Landau level 
at the charge neutrality point is 8 -fold degenerate con- 
sisting of 4 electron and 4 hole levels. The other Landau 
levels are each 4-fold degenerate and have a hole ("-" 
sign) or electron ("+" sign) character. The density of 
states of each Landau level consists of a region of ex- 
tended states, which carry current, near the center, while 
states in the upper or lower tails are localised. Minima 
of p xx occur when the Fermi energy (which can be var- 
ied by sweeping gate voltage or field) lies in the middle 
between two Landau levels. The minima, correspond- 
ing to v = ±4 and v = +8, are shown in Fig. |2jb) by 
straight arrows. 

To determine activation gaps we measured the lon- 
gitudinal resistivity p xx as a function of charge car- 
riers concentration n at various temperatures T - 
0.5 K-220 K Fig. 13b). At the highest field used, Shub- 
nikov de Haas oscillations remain visible up to 220 K. 
The resistivity at the minima remains zero up to 24 K. 
This temperature is considerably low compared with 
single layer graphene, where the quantum Hall effect 
can even be observed up to room temperature II 1311 . 

With increasing temperature electrons can be ther- 
mally activated from localised states to the extended 
states of the neighboring Landau levels (arrows in the 
Fig. [3J a )) increasing p xx at a minimum. The temper- 
ature dependences of p xx for v = 4 and v = 8 are 
shown in Fig. |2jc) and (d) (at v = -4 and v = — 8 
similar results are obtained). For low temperatures the 
dependence of p xx on temperature flattens off signifi- 
cantly, which is commonly attributed to non-activated 
transport mechanism as variable range hopping |Qi" 
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Figure 2: Density of states modelled with extended (black areas) and 
localised (dashed areas) states (a). Longitudinal resistivity p xx as a 
function of charge carrier concentration at B = 15 T and various tem- 
peratures (b). Temperature dependences of p xx at various fields at 
v = 4 (c) and v = 8 (d). The solid lines represent a fit, using the 
Fermi-Dirac distribution function. 



The gray lines in Fig.|2jc), (d) show a fit of the exper- 
imental p xx (T) dependences by Fermi-Dirac function 
1/(1 + exp(A a /2lcBT)) to describe the activated trans- 
port. The activation gaps A fl are used as a fit parameter. 

We will show below that the activation gaps measured 
do not exceed the thermal energy considerably. There- 
fore we use the Fermi-Dirac function to extract the gaps, 
rather than the commonly used low temperature expo- 
nential approximation (Arrhenius plot). We assume that 
the conductivity is determined by the thermal activa- 
tion to the single energy between localised and extended 
states, namely the conductivity edge Eq of the Landau 
levels (Fig. [2ja)). The thermal activation deeper inside 
the Landau levels is neglected. We have verified the va- 
lidity of this assumption by modelling the longitudinal 
resistivity p xx inside a gap A fl between two broadened 
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Landau levels using standard Kubo-Greenwood formal- 
ism 1 15, T^. Indeed, for A a > 2k[,T the gaps given by 
this model agree within 2% with those extracted from 
the Fermi-Dirac fit. 

The experimental gaps, extracted from the Fermi- 
Dirac fit are presented in Fig. [3] Due to the flattening of 
experimental p xx {T) curves at low temperatures (T < 70 
K for v = +4 and T < 50 K for v = +8) we limited our 
fitting range to higher temperatures but included a non- 
activated small constant background. Extending the fit- 
ting to lower temperature yields systematically different 
gaps. We indicated this effect by the size of the error 
bars in Fig. [3] 

The activation gap determined this way corresponds 
to the double energy from the Fermi energy to the near- 
est extended state and is therefore given by A„ = AE-T, 
with a Landau level broadening F and an energy AE be- 
tween the centers of the levels . The measured energy is 
thus lower than the bare gap of infinitely sharp Landau 
levels and even vanishes when F = AE. 

Theoretically the low band Landau levels in bilayer 
graphene are given by llT 
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with N = 0, ±1, ±2 . . ., the "+" sign corresponds to 
electrons or holes. In this equation j\ = 0.4 eV lfl8l[T9 1 
is the interlayer coupling constant and c — 10 6 ms~' fl9[] 
is the velocity in single layer graphene. For low mag- 
netic field Eq. (Q} reduces to yfl: 
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where the effective mass arises from the interlayer cou- 
pling j\ as m* = j\l(2c 2 ). The low field expression 
Eq. (f2]i predicts a linear field dependence, arising from 
a parabolic dispersion, with different slopes depending 
on N. The full Eq. (Q]) also displays a non-parabolicity 
for high fields and/or energies. 

In Fig. [3] we compare our experimental data with the 
calculated energy splitting for v = ±4 and v = ±8 using 
Eqs. ([TJ and (0. From the figure it can be seen that for 
the same magnetic field the experimental gaps for the 
v = +8 (N — 1 — > 2) are always smaller than that for 
v = ±4 (N = — > 1) as theoretically predicted. Fur- 
thermore the activation gaps seem to have a sublinear 
field dependence which is more pronounced for v = ±8 
and which is partially attributed to the predicted non- 
parabolicity (Eq. ((TJ)- The experimental gaps vanish at 
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Figure 3: (color online). The activation gaps A a extracted from the fit 
(Fig-EJc), (d)) as a function of magnetic field at v = ±4 (a) and v = ±8 
(b). Black filled (blue empty) circles correspond to electrons (holes). 
The black solid lines show the distances between centers of corre- 
sponding Landau levels as given by hyperbolic dispersion Eq. Q, the 
gray solid lines represent the parabolic approximation Eq.{2). The 
dashed lines represent the expected field dependence for a constant 
level broadening, the dash-dotted line includes a field dependent level 
width (see the main text for details). 



a finite field below 10 T which we attribute to the van- 
ishing of the activation gaps when broadened Landau 
levels start to overlap, i. e. when T = AE. For that rea- 
son, the measured activation gaps follow the bare gaps 
with increasing field with an offset that is determined by 
r. This r appears to be a constant To for the v = ±4 
transition (dashed line in Fig. |3ja)) whereas it shows 
a field dependence F = To + aB for v = +8 transi- 
tion (dashed-dotted and dashed lines in Fig. Ob)). This 
field dependent broadening for the higher Landau lev- 
els (N + 0) could be caused by the rippled nature of 
the graphene, which leads to a difference of the perpen- 
dicular component of the magnetic field throughout the 
sample and therefore to an inhomogeneous linearly de- 
pendent broadening. This field dependent mechanism 
was also observed for non-zero energy Landau level in 
single layer graphene B20I1 . The lowest Landau level 
(N = 0) is not field dependent and thus the N = — > 1 
transition is much less affected by this mechanism. 

It is interesting to see, that for bilayer graphene our 
measurements show a broadened zeroth Landau level up 
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to 30 T, without any indication of a narrowing as was 
seen in single layer graphene I20I1 . It is possible that 
we need samples of a better quality or higher magnetic 
fields to observe this narrowing in bilayer graphene. It 
might also be an indication of a different broadening 
mechanism of the lowest Landau level in the single and 
the bilayer graphene and therefore different conditions 
for the coexistence of holes and electrons at the charge 
neutrality point. 

3. Conclusions 

To conclude, we have investigated the gaps between 
Landau levels at filling factors v = ±4 and v = +8 in 
bilayer graphene by means of activated transport mea- 
surements. Experimental field dependences follow Lan- 
dau level quantization for a 4-band hyperbolic relativis- 
tic dispersion with field dependent offsets caused by a 
field dependent level width. 
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